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We propose a nonextensive generalization (q parametrized) 
of the von Neumann equation for the density operator. Our 
model naturally leads to the phenomenon of decoherence, and 
unitary evolution is recovered in the limit of q — > 1. The 
resulting evolution yields a nonexponential decay for quantum 
coherences, fact that might be attributed to nonextensivity. 
We discuss, as an example, the loss of coherence observed in 
trapped ions. 

05.30.Ch, 03.65. Bz, 42.50.Lc 



In the past decade, there have been substantial ad- 
vances regarding a nonextensive, g-parametrized gener- 
alization of Gibbs-Boltzmann statistical mechanics . 
The q parametrization is based on an approximate ex- 
pression for the exponential, or the g-exponential 



e q (x) = [1 + (1 - «): 



l V(l-9) 



(1) 



being the ordinary exponential recovered in the limit 
of q — ► 1 (ei(ar) = e x ). Thus we have a nonex- 
tensive g-exponential, i.e., e q (x)e q (y) ^ e q (x + y), in e -g 
general. Such a parametrization is the basis of Tsal- 
lis statistics, in which it is used a g-parametrized nat- 
ural logarithm for the definition of a generalized en- 
tropy. We recall that in Tsallis' statistics the entropy 
follows the pseudo- additive rule |l],^| S q (A + B) = S q (A) 
+ S q {B) + (1 - q)S q (A)S q (B), where A and B represent 
two independent systems, i.e., Pij(A + B) — pi(A)pj(B). 
Hence for q < 1 we have the supper-additive regime, and 
for q > 1 the sub-additive regime. The parameter q thus 
characterizes the degree of extensivity of the physical sys- 
tem considered. Tsallis' entropy seems to be previligied, 
in the sense that it satisfies the concavity criterion, i.e., 
it has a well defined concavity Moreover, it has 

also been shown that most of the formal structure of 
the standard statistical mechanics (and thermodynam- 
ics) is retained within the nonextensive statistics, e.g., 
the H-theorem IjJ, and the Onsager reciprocity theorem 
[^J. Nonextensive effects are usual in many systems ex- 
hibiting, for instance, long-range interactions, and the 
nonextensive formalism has been succesfully applied not 
only to several interesting physical problems, but also 
outside physics 0]. We may cite for instance, applica- 
tions to statistical mechanics [||, Levy anomalous sup- 
perdiffusion J?j], quantum scattering pj, low dimensional 
dissipative systems |J , in cosmology , and others [^) . 
This means that departures from exponential (or loga- 
rithmic) behavior not so rarely occur in nature, and that 



the parametrization given by Tsallis' statistics seems to 
be adequate for treating them. Discussions concerning 
the quantification of quantum entanglement pd| , ^2[ as 
well as the implications on local realism p"3| ] may be 
also found, within the nonextensive formalism. Quan- 
tum entanglement has itself a nonlocal nature, and this 
may have somehow a connection with the general idea of 
nonextensivity. In fact, it has been concluded that entan- 
glement may be enhanced in the supper-additive regime 
(for q < 1) |pd| , ^2| . Thus, although nonextensive ideas 
open up new interesting possibilities, their application 
to a field such as foundations of quantum mechanics has 
not been so frequent. We would like, in this letter, to 
address fundamental aspects of quantum theory within a 
nonextensive approach. 

In quantum mechanics, the issue of decoherence, or the 
transformation of pure states into statistical mixtures has 
been atractting a lot attention, recently, especially due 
to the potential applications, such as quantum comput- 
ing and quantum criptography j^J, that might arise in 
highly controlled purely quantum mechanical systems, 
e.g., in trapped ion systems p5[ . That extraordinary con- 
trol is also allowing to address fundamental problems in 
quantum mechanics, such as the Schrodinger cat prob- 
lem 16 1, and the question of the origin of decoherence 
itself. Despite of the progresses, it does not exist yet a 
proper theory handling the question of the loss of coher- 
ence in quantum mechanics, although there are several 
propositions, involving either some coupling with an en- 
vironment [|l7 18 1, or spontaneous (intrinsic) mechanisms 
ifLsf . Indeed dissipative environments, in which it is veri- 
fied loss of energy, tend to cause decoherence. In spite of 
the destructive action of dissipation, though, it has been 
worked out a scheme to recover quantum coherences in 
cavity fields JljJ, even if its environment is at T ^ 0, 
and after the system's energy, and coherences, have sub- 
stantially decayed. Regarding the intrinsic decoherence, 
several models have been already presented Jl8| , and de- 
coherence has been attributed, for instance, to stochastic 
jumps in the wave function |pof , or gravitational effects 
|^lf . These models may contain one |po| or even two new 
parameters It is normally proposed some kind of 

modification of the von Neumann equation for the den- 
sity operator p 



dp 



[H,p]- 



(2) 



A typical model yielding intrinsic decoherence, such as 
Milburn's p(| (see also reference p3] ) , gives the following 



1 



modified equation for the evolution of p 



dp 
~dt 



J[^/5]-^2^, [H,p]}, 



(3) 



where r is a fundamental time step. The second term in 
the right-hand-side (double commutator), leads to an ex- 
ponential decay of coherences without energy loss |^0|, 24 1 . 
Nevertheless, other models concerning decoherence [21|, 
as well as quantum measurements p5fl , predict nonexpo- 
nential decay of coherences, going as exp(— jt 2 ) rather 
than as exp(— jt). In plj , decoherence is attributed to 
a whormhole-matter coupling of nonlocal nature. This 
suggests that decoherence might not be appropriately 
described by a Markovian stochastic process (2^] . More- 
over, experimental evidence of nonexponential decay of 
quantum coherences has been recently reported, in a ex- 
periments involving trapped ions, J2(|, as well as in quan- 
tum tunneling of cold atoms in an accelerating optical 
potential p7fl . 

Here we propose a novel model to treat decoherence, 
within a nonextensive realm. We are able to write down 
a g-parametrized (generalized) evolution for the density 
operator, which naturally leads to decoherence. As a first 
step we may express von Neumann's equation [Eq. (j^)] 
in terms of the superoperator C, being ihCp = Hp — pH, 
so that we may write its formal solution in a simple form 



pit) = cMct)p(o). 



(4) 



Now, we substitute the exponential superoperator 
exp(£t) by a q-parametrized exponential as defined in 
Eq. (§ 



- 1 

p(t)= l + (l-q)a p(0). 



(5) 



The standard unitary evolution given by von Neumann's 
equation is recovered in the limit of q — > 1. Hence, de- 
parting from Eq. (|^), we may write a generalized (q- 
paramctrized) von Neumann equation, or 



dp 
~dl 



C 



I + (1 - q)Lt 



Pit)- 



(6) 



Here we consider the case in which the parameter q is 
very close to one. This represents a small deviation from 
the unitary evolution given by (|^), but yields loss of co- 
herence. For |1 — q\flt <C 1, where fi is the characteris- 
tic frequency of the physical system considered, we may 
write 



dp 
~dl 



l-{l-q)tl 2 p(t). 



(7) 



It may be shown, from Eq. (J?]), that there is decay of the 
nondiagonal elements in the basis formed by the eigen- 
states of H, characterizing decoherence, while the diag- 
onal elements are not modified. Therefore the density 



operator remains normalized at all times, or Trp(t) = 
Trp(0) = 1, property which must hold in any proper 
model for decoherence. We shall remark that in order to 
have a physically acceptable dynamics, the extensivity 
parameter should be greater than one (q > 1), i.e., in the 
sub-additive regime. We may also rewrite Eq. (fij) as 



dp 
~dt 



[H,p]+g(t)[H, [H,p]}, 



(8) 



where g(t) = (1 - q)t/h 2 . Eq. (g) is similar to Eq. (g), 
apart from the time-dependent factor g(t) multiplying 
the double commutator. Note that Eq. (g) is only valid 
for short times, or |1 — q\flt <C 1. We have therefore 
constructed a novel model, based on a nonextensive gen- 
eralization of von Neumann's equation, which predicts 
loss of coherence, although with a time-dependence for 
the nondiagonal elements in the energy basis different 
from most that can be found in the literature |18|. This 
particular time-dependence will lead to a nonexponential 
decay of quantum coherences, as we are going to discuss. 

Now we would like to apply our model to a specific sys- 
tem in which experimental results are available. One of 
the most interesting systems investigated nowadays is a 
single trapped ion interacting with laser fields. Quantum 
state engineering of motional states of a massive oscil- 
lator has been achieved PH , and the loss of coherence 
in that system has been already observed |2^,|26|]. In 
the experiment described in reference pq j, two Raman 
lasers induce transitions between two internal electronic 
states (|e) and \g)) of a single 9 Be + trapped ion, as well 
as amongst center-of-mass vibrational states \n). In the 
interaction picture, and under the rotating wave approxi- 
mation, the effective hamiltonian (in one dimension) may 
be written as lEsi 



H 



being f2 the coupling constant, 6 the detuning of the fre- 
quency difference of the two laser beams with respect 
to wo, which is the energy difference between |e), or 
\g) (E e — Eg = Tiluq). The Lamb-Dicke parameter is 
•q — kyfh/ (2moj x ), where k is the magnitude of the dif- 
ference between the wavevectors of the two laser beams, 
lo x is the vibrational frequency of the ion's center of mass 
having a mass m, and d x {d\.), are annihilation and cre- 
ation operators of vibrational excitations. The ion is Ra- 
man cooled to the (vibrational) vacuum state |0), and 
from that Fock (as well as coherent and squeezed) states 
may be prepared by applying a sequence of laser pulses 
p8[ . If the atom is initially in the \g, n) state, and the Ra- 
man lasers are tuned to the first blue sideband (S — w x ), 
the evolution according to the hamiltonian in Eq. @ (for 
small if) will be such that the system will perform Rabi 
oscillations between the states | g, n) and |e, n + 1). If the 
excitation distribution of the initial vibrational state of 



2 



the ion is P n , occupation probability of the ground state 
P g (t) will be given by 



P g (n,t) = 



1 



1 + ^2 P n cos(2fi„t) e" 



Here the Rabi frequency Q n is 



(10) 



(11) 



where Q,/2tt = 500kHz, are Laguerre generalized poly- 
nomials, and 7„ = 7o(n + l) ' 7 , with 70 = 11.9kHz. 

The experimental results are fitted using Eq. (|H 
pj| . The Rabi oscillations are damped, and an empiri- 
cal damping factor, 7„, is introduced in order to fit the 
experimental data Q . There heve been attempts to de- 
rive such an unusual dependence on n for the damping 
factor 7„, by taking into account the fluctuations in laser 
fields and the trap parameters j2j|. Those models are 
somehow connected to a evolution of the type given by 
Eq. (|3|) |2(],Q , which predicts loss of coherence without 
loss of energy. In fact, at a time-scale in which the ion's 
energy remains almost constant, decoherence is consider- 
ably large Jl5|,|28| . Moreover, the actual causes of loss of 
coherence in experiments with trapped ions are still not 
identified |]l5| , ^o| . Our approach differs from previous 
ones because we depart from a different dynamics which 
leads to a peculiar nonexponential time-dependence. By 
employing a similar methodology to the one discussed in 
reference , we may calculate, from Eq. (^) , the prob- 
ability of the atom to occupy the ground state, obtaining 



cos(2f2 ra t) e 



(12) 



where j niq = (q — l)fi„/2, i.e., the damping factor aris- 
ing in our model depends on the Rabi frequency fl n and 
on the parameter q in a simple way. It is also verified 
a nonexponential decay which goes as exp(— jt 2 ) rather 
than an exponential one. Now we may proceed with a 
graphical comparison between the curves plotted from 
expressions (|l0|) and (|l2|). This is shown in Figure 1, for 
two distinct cases. In Figure la we have plotted both 
the probability P g in Eq. (jlj) (dashed line) and the 
one in Eq. ( |lfj| ) (full line), as a function of time, for an 
initial state \g,Q) (P n — ^n,o) 5 a Lamb-Dicke parame- 
ter r] = 0.202, and q = 1.001. We note a clear depar- 
ture from exponential behavior in the curve arising from 
our model, although they might coincide for some time- 
intervals. We remark that Eq. (12) is valid for times 
such that \l-q\flt < 1. Here \ l-q\flt max « 0.17. Right 
below, in Figure lb, there is a similar plot, but using 
a different initial condition for the distribution of exci- 
tations of the ion state, e.g., a coherent state \a), with 
n = a 2 — 3 (P n — exp(— a 2 )a 2n /n\) instead of the vac- 
uum state |0). In the coherent state case both curves 



are very close, although it seems that our curve (dashed 
line) represents a better fit for the experimental results 
p8| |. In general, our results are in reasonable agreement 
with the experimental data, as it may be seen in Figure 1 
and in reference | p8| . This means that a nonexponential 
decay may also account for the loss of coherence in the 
trapped ions system. 
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FIG. 1. Probability of having the atom in the ground state 
P g as a function of time; (a) for an initial (vibrational) vacuum 
state and (b) for an initial (vibrational) coherent state with 
n = 3. In either case rj = 0.202, and n/2n = 500kHz. The 
solid lines are calculated from from Eq.(10), and the dashed 
lines from Eq.(12), with an extensivity parameter q — 1.001. 

In our model, the effects leading to decoherence are 
embodied within the nonextensive parameter q, rather 
than into a "fundamental time step" r, as it occurs in 
Milburn's model. Nonextensivity may be especially rel- 
evant when nonlocal effects, for instance, long-range in- 
teractions, or memory effects are involved [||. In the 
specific example treated here, the ion is not completely 
isolated from its surroundings; electric fields generated in 
the trap electrodes couple to the ion's charge, and this is 
considered to be a genuine source of decoherence to the 
vibrational motion of the ion Jll||3(J . 

In conclusion, we have presented a novel approach for 
treating the loss of coherence in quantum mechanics, 
based on a nonextensive formalism. In our model, the 
quality of decoherence depends on a single parameter, q, 
related to the nonextensive properties of the physical sys- 
tem considered. We obtain, with such a parametrization, 
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an evolution equation which is a (nonextensive) general- 
ization of von Neumann's equation for the density opera- 
tor, and that leads, in general, to a nonexponential decay 
of quantum coherences. We have applied our model to 
a concrete physical problem, that is, the decoherence oc- 
curring in the quantum state of a single trapped ion un- 
dergoing harmonic oscillations. This is a well know case 
in which decoherence may be readily tracked down, yet 
its causes remain unexplained. We have found that for 
values of the parameter q rather close to one (in the sub- 
additive regime), our model is in reasonable agreement 
with the available experimental results, without the need 
of introducing any supplementary assumptions. 
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